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Abstract 

We calculate spectral functions associated with hadronic current correlation functions for vector 
currents at finite temperature. We make use of a model with chiral symmetry, temperature- 
dependent coupling constants and temperature-dependent momentum cutoff parameters. Our 
model has two parameters which are used to fix the magnitude and position of the large peak seen 
in the spectral functions. In our earlier work, good fits were obtained for the spectral functions that 
were extracted from lattice data by means of the maximum entropy method (MEM). In the present 
work we extend our calculations and provide values for the three-momentum dependence of the 
vector correlation function at T = 1.5Tc. These results are used to obtain the correlation function 
in coordinate space, which is usually parametrized in terms of a screening mass. Our results for 
the three-momentum dependence of the spectral functions are similar to those found in a recent 
lattice QCD calculation for charmonium [S. Datta, F. Karsch, P. Petreczky and I. Wetzorke, 
hep-lat/0312037 . However, we do not find the expontential behavior in coordinate space that 
is usually assumed for the spectral function for T > Tc and which allows for the definition of a 
screening mass. 
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I. INTRODUCTION 



In a number of recent works [ 1- 3] we have calculated various hadronic correlation func- 
tions and compared our results to results obtained in lattice simulations of QCD [ 4- 6] . The 
lattice results for the correlators, G{t, T), may be used to obtain the corresponding spectral 
functions, a{u!,T), by making use of the relation 

/•oo 

G{t,T)^ / duja{uj,T)K{T,uj,T), (1.1) 
Jo 

where 

^ ' ' ^ sinh(cc;/2T) ^ ^ 

The procedure to obtain a{uj, T) from the knowledge of G(t, T) makes use of the maximum 
entropy method (MEM) [7-9], since G{t,T) is only known at a hmited number of points. 

In our studies of meson spectra at T = and at T < Tc we have made use of the Nambu- 
Jona-Lasinio (NJL) model. The Lagrangian of the generalized NJL model we have used in 
our studies is 



C = q{ij - m')q + ^ Y^^qX'qf + {qil^^qf] (1-3) 

i=0 

G 

+ ^{det[g(l + As)?] + det[g(l - As)^]} + C^onj • 

Here, mP is a current quark mass matrix, mP — diag{m^,m'^,m'^). The Aj are the Gell- 
Mann (flavor) matrices and A° = \/2/31, with 1 being the unit matrix. The fourth term is 
the 't Hooft interaction and Cconf represents the model of confinement used in our studies 
of meson properties. 

In the study of hadronic current correlators it is important to use a model which respects 
chiral symmetry, when mP — 0. Therefore, we make use of the Lagrangian of Eq. (1.3), 
while neglecting the 't Hooft interaction and Cconf - In order to make contact with the results 
of lattice simulations we use the model with the number of flavors, Nf ~ 1. Therefore, the 
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A* matrices in Eq. (1.3) may be replaced by unity. We then use 

L = q{ii - m')q + '^[{qqf + [qil.qf] (1.4) 

in order to calculate the hadronic current correlation functions. Thus, there are essentially 
three parameters to consider, Gs, Gy and a Gaussian cutoff parameter a, which restricts 
the momentum integrals through a factor exp[— A; ^/«^]. When we use the NJL model 
to study matter at finite temperature, we introduce the temperature-dependent parameters 
Gs{T), Gv{T) and a{T). These parameters have been adjusted to obtain fits to the spectral 
functions a{ijj,T) for T /Tc — 1.5 and T/Tc — 3.0, which are the values of T/Tc studied in 
the lattice simulations of QCD [10]. 

The temperature-dependent coupling constants and cutoff parameters of our work are 
analogous to the corresponding density- dependent parameters introduced in Ref. [11] and 
[12]. Further study of models with temperature-dependent and density-dependent parame- 
ters are of interest and a general theoretical formalism for the introduction of such depen- 
dencies should be considered. 

In this work we limit our study to the data for the vector correlator at T = l.STc [10, 
13] and therefore only need to specify Gy(T) and oi{T) at that temperature. In Figs. 1 
and 2 we show the data obtained by the MEM method at T/Tc = 1-5 and T/Tc = 3.0 for 
both pseudoscalar and vector correlators [10, 13]. The second peaks in these correlators are 
known to be a lattice artifacts [13] and our calculations do not reproduce those peaks. 

The organization of our work is as follows. In Appendix A we present the formalism for 
calculation of pseudoscalar and vector correlation functions making use of the Lagrangian of 
Eq. (1.4). In Appendix A we discuss the calculation of the correlator in the case the quark 
and antiquark carry zero total momentum. In Appendix B we show how the formalism is 
modified for the correlator calculated at finite momentum P. In Section II we present the 
results of our calculation of the imaginary part of the correlator a{uj, P). (Since we place P 
along the z-&xSs this quantity may be written as (7(0;, 0, 0, P^ in accord with the notation of 
Ref. [10].) Our results for a{uj, P) will be presented for a series of values of | P| in Section 
II. In Section III we present our result for the coordinate-dependent correlator C{z) which 
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is proportional to the correlator defined in Eq. (1) of Ref. [10], 

cw=ir.P.e-rd.^(^M^. (1.6) 

^ J-oo Jo ^ 

We may also use the form 

cw^i£.py-/^P^!:(fl^. (1.6) 

Finally, in Section IV we present our conclusion and further discussion. 

II. TEMPERATURE DEPENDENT HADRONIC CURRENT CORRELATORS 
AT FINITE MOMENTUM 

We make use of the formalism presented in Appendices A and B to obtain values of the 
vector correlator at T = l.bTc- Here we take Tc — 270 MeV, since we are considering lattice 
calculations made in the quenched approximation. 

In Fig. 3 we present (y{uj)/u'^, for various values of |P|, as function of uj/T. Comparison 
may be made to the lattice data shown in Fig. 2 [10, 13]. (We again note that our calculation 
does not reproduce the second peak in the lattice data which is known to be a lattice artifact 
[13].) The curves shown in Fig. 3 are given for values of | P | ranging from 0.10 GeV to 3.10 
GeV in steps of 0.20 GeV. The same data are presented in Fig. 4, where our results for 
a{(jj)/uj'^ are given as a function of a;^. In this figure we see that the maximum value of o;^ 
used in our calculations is 90 GeV^, so that ujmax = 9.49 GeV. Our results for the various 
values of |P| given in Fig. 3 may be compared to Fig. 20 of Ref. [14]. We see that the 
results calculated by completely different methods are similar. 



III. THE COORDINATE-SPACE CORRELATION FUNCTION 

We have used the results of our calculations which were presented in Fig. 4 to calculate 
C{z) of Eq. (1.6). The result of that calculation is shown in Fig. 5. We note that the simple 
assumption for the behavior of this correlator that is usually made, C{z) ~ exp[— msc-z], is 
not born out in our calculation. Therefore, we do not provide a value for the screening mass 
rUsc which is defined when assuming exponential behavior of the correlator. For the study 
of charmonium on the lattice, values found for the screening mass are given in Ref. [14]. 
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FIG. 1: The spectral functions cr/o;^ for pseudoscalar states obtained by MEM are shown [10, 13]. 

The sohd Unc is for T /T^ = 1.5 and tlic daslicd Unc is for T/Tc = 3.0. The second peak is tliought 
to be a lattice artifact [13] and is not reproduced in our model. (See Fig. 3, for example.) 

0.18 -, 




co/T 

FIG. 2: The spectral functions a/uP' for vector states obtained by MEM are shown [10, 13]. See 
the caption of Fig. 1. These results may be compared to the theoretical results shown in Fig. 3 
where we see the absence of the second peak. That peak is thought to be a lattice artifact [13]. 
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0.18-, 




co/T 

FIG. 3: The imaginary part of the correlator a{io)luP' is shown for various values of \P\ as a 
function of lo/T. Starting with the topmost curve the values of \P\ in GeV units are 0.10, 0.30, 
0.50, 0.70, 0.90, 1.10, 1.30, 1.50, 1.70, 1.90, 2.10, 2.30, 2.50, 2.70, 2.90, and 3.10. Here we have 
used Gs = 0.66 GeV-^ and a = 4.4 GeV. 

0.18-1 




20 40 60 80 100 

(0 (GeV') 

FIG. 4: The imaginary part of the vector correlation function a{u:)/oj'^ is shown for various values 
of |P| as a function of lo"^. See the caption of Fig. 3. 
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FIG. 5: The correlation function C{z) defined in Eq. (1.6) is shown. 
IV. DISCUSSION 

Recent theoretical work concerning the quark-gluon plasma has been discussed by 
Shuryak [15]. He notes that the physics of excited matter produced in heayy ion colli- 
sion in the region T^. < T < AT^ is different from that of a weakly coupled quark-gluon 
plasma. That is due to the strong coupling generated by the bound states of the quasipar- 
ticles. These bound states appear as resonant structures when the imaginary parts of the 
hadronic correlators are extracted from lattice data using the maximum entropy method 
(MEM) [4-9] . These resonances lead to very strong interactions between the quasiparticles 
and are, in part, responsible for very small mean free paths and collective flow in heavy- 
ion collisions. That flow may be described by hydrodynamics. Indeed, Shuryak suggests 
that "... if the system is macroscopically large, then its description via thermodynamics of 
its bulk properties (like matter composition) and hydrodynamics for space-time evolution 
should work" [15]. 

In our present work we have presented a simple chiral model that is able to reproduce the 
resonances extracted from the lattice data via the MEM procedure. We have calculated the 
imaginary parts of the vector correlator for various values of the total external momentum 

— * 

P. We have also calculated the correlation function, C{z), and do not find the simple 
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exponential behavior C{z) ~ exp[— msc-z] usually assumed for that correlator. In the special 
case in which the correlator is dominated by a simple delta function for small cu, we may 
define m^c and see that the screening and pole masses are equal [10].) 

APPENDIX A 

For ease of reference, we present a discussion of our calculation of hadronic current 
correlators taken from Ref . [ 3 ] . The procedure we adopt is based upon the real-time finite- 
temperature formalism, in which the imaginary part of the polarization function may be 
calculated. Then, the real part of the function is obtained using a dispersion relation. The 
result we need for this work has been already given in the work of Kobes and Semenoff [ 16 ] . 
(In Ref. [ 16 ] the quark momentum is k and the antiquark momentum is k — P. We will adopt 
that notation in this section for ease of reference to the results presented in Ref. [ 16 ].) With 
reference to Eq. (5.4) of Ref. [16], we write the imaginary part of the scalar polarization 
function as 

ImJsiP^T) = -N,Pse{P') /-^e-*^V«^ ( r^^-jTr] (Al) 

x{[l - ni(fc) - n2(fc)](5(P° - Er{k) - E^ik)) 
-[ni(fc) - n2(A;)]5(P° + E^{k) - E^^k)) 
-[n2{k) - m{k)]5{P'' - Ei{k) + E2{k)) 
-[1 - n^{k) - n2(A;)]5(P° + E^{k) + E^ik))} . 

Here, Ei{k) = [P + ml{T) ] V^. Relative to Eq. (5.4) of Ref. [ 16 ], we have changed the sign, 
removed a factor of g'^ and have included a statistical factor of Nc- In addition, we have 
included a Gaussian regulator, exp[ — A; ^/a^ ]. The value a = 0.605 GeV was used in our 
applications of the NJL model in the calculation of meson properties at T = 0. We also 
note that 

and 
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For the calculation of the imaginary part of the polarization function, we may put — 
ml{T) and {k — Py — ml{T), since in that calculation the quark and antiquark are on- 
mass-shell. In Eq. (Al) the factor Ps arises from a trace involving Dirac matrices, such 
that 

(3s = -Tr[(^ + mi)(^-f + m2)] (A4) 
= _2(mi + m2)^ (A5) 

— * 

where mi and 1712 depend upon temperature. In the frame where P = 0, and in the case 
mi = m2, we have (3s = 2Pq^(1 — Am'^ / Pq). For the scalar case, with mi = m2, we find 

?2 / A^2frr\ \ 3/2 



,..,(P^^) = ^(l-^!^P)'e-/-[:-2M*)l. 



(A6) 



where 



p ^^l-m'iT). (A7) 
For pseudoscalar mesons, we replace (3s by 

(3p = -'B:[i75(|^ + mi)i75(|^ - f + ^2) ] (A8) 

= - 2(mi - m2)^ , (A9) 

which for mi = m2 is (3p — 2Pq in the frame where P = 0. We find, for the tt mesons, 

Im Jp(P^ T) = ^ (1 - e-^^ [ 1 _ 2n,{k) ] , (AlO) 

where k"^ — Pq/4 — m^{T), as above. Thus, we see that the phase space factor has an 
exponent of 1/2 corresponding to a s-wave amplitude. For the scalars, the exponent of the 
phase-space factor is 3/2, as seen in Eq. (A6). 
For a study of vector mesons we consider 

P^. = + - f + m2) ] , (All) 

and calculate 

gt'^P^^ =4[P^ -ml-ml + 4mim2 ] , (A12) 

which, in the equal- mass case, is equal to 4Po -|- 8m^(T), when mi = m2 and P — 0. This 
result will be needed when we calculate the correlator of vector currents. Note that, for 
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the elevated temperatures considered in this work, mu{T) — md{T) is quite small, so that 
4Po^ + 8m^(T) can be approximated by 4P|, when we consider the vector current correlation 
functions. In that case, we have 

Im (P^ T) ~ ^Im Jp{P\ T) . (A13) 



At this point it is useful to define functions that do not contain that Gaussian regulator: 

Im Jp(P^ T) = ^ fl - 2n,{k) ] , (A14) 







and 



Im Jy{P\ T) = ( 1 - 2n,{k) ] , (A15) 



3 Stt V' Pq / 

For the functions defined in Eq. (A14) and (A15) we need to use a twice-subtracted dispersion 
relation to obtain Re Jp{P^, T), or Re JviP"^, T). For example. 

Re Jp(P^ T) = Re Jp(0, T) + [ Re Jp(Po , T) - Re Jp(0, T) ] (A16) 



, P\P' - Pj) r Im Jp(s,T) 



where A^ can be quite large, since the integral over the imaginary part of the polarization 
function is now convergent. We may introduce Jp{P'^,T) and Jv{P^,T) as complex func- 
tions, since we now have both the real and imaginary parts of these functions. We note that 
the construction of either Re Jp(P^,T), or Re Ji/(P^,T), by means of a dispersion relation 
does not require a subtraction. We use these functions to define the complex functions 
Jp{P\T) and Jv(P^T). 

In order to make use of Eq. (A16), we need to specify Jp(0) and Jp{Pq). We found it 
useful to take P| = -1.0 GeV^ and to put Jp(0) = Jp(0) and Jp(Po^) = Jp{Pq)- The 
quantities Jy(0) and Jv{Pq) are determined in an analogous function. This procedure in 
which we fix the behavior of a function such as ReJy(P^) or ReJy(P^) is quite analogous 
to the procedure used in Ref. [17]. In that work wc made use of dispersion relations to 
construct a continuous vector-isovector current correlation function which had the correct 
perturbative behavior for large P^ ^ — oo and also described the low-energy resonance 
present in the correlator due to the excitation of the p meson. In Ref. [ 17 ] the NJL model 
was shown to provide a quite satisfactory description of the low-energy resonant behavior 
of the vector-isovector correlation function. 
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We now consider the calculation of temperature-dependent hadronic current correlation 
functions. The general form of the correlator is a transform of a time-ordered product of 
currents, 

^C(P^ T)^J d'xe'''-^ « T{j{x)m) » , (A17) 

where the double bracket is a reminder that we are considering the finite temperature case. 

For the study of pseudoscalar states, we may consider currents of the form jp,i{x) — 
q{x)i'-f5X^q{x), where, in the case of the tt mesons, i = 1,2 and 3. For the study of scalar- 
isoscalar mesons, wc introduce js,i{x) = q{x)X'^q{x), where i = for the flavor-singlet current 
and i = 8 for the flavor-octet current [18]. 

In the case of the pseudoscalar-isovector mesons, the correlator may be expressed in terms 
of the basic vacuum polarization function of the NJL model, Jp{P'^,T) [18-20]. Thus, 

CAP\T)^MP\T)^-^^-^±j-^, (A18) 

where Gp{T) is the coupling constant appropriate for our study of tt mesons. We have 
found Gp{T) — 13.49 GeV~^ by fitting the pion mass in a calculation made at T = 0, with 
f^u — iTT'd — 0.364 GeV. The result given in Eq. (A18) is only expected to be useful for small 
P^, since the Gaussian regulator strongly modifies the large behavior. Therefore, we 
suggest that the following form is useful, if we are to consider the larger values of P^. 



CpjP^T) 

p2 



Jp{P\T) 

p2 



^ (A19) 



l-Gp{T)Jp{P^,T) 



(As usual, we put P = 0.) This form has two important features. At large Pq, 
lmCp(Po, T)/Pq^ is a constant, since ImJp(Pg^,T) is proportional to P^. Further, the 
denominator of Eq. (A19) goes to 1 for large Pq . On the other hand, at small Pq , the 
denominator is capable of describing resonant enhancement of the correlation function. As 
we have seen, the results obtained when Eq. (A19) is used appear quite satisfactory. (We 
may again refer to Ref. [17], in which a similar approximation is described.) 

For a study of the vector-isovector correlators, we introduce conserved vector currents 
~ Q{^)'liJ'^iQ{^) with i=l, 2 and 3. In this case we define 

J^(P^T) = (^^^'^ - Jv{P\T) (A20) 
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and 



(A21) 



taking into account the fact that the current j,j,,i{x) is conserved. We may then use the fact 
that 

1 



and 



Im Jv{P^T) 



P2 + 2ml{T) 



Stt 



-ImJp(P^T). 



AmliT) 



p2 



1/2 



(See Eq. (A7) for the specification of A; = |A;|.) We then have 

Cv{P',T) = MP',T)— ^ 
where we have introduced 



Gv{T)MP^T)' 



(A22) 

e-W[i_2ni(A;)] (A23) 

(A24) 

(A25) 



lmJv{P\T) 



Stt 



^ImJp(P^T). 



(l-— [l-2ni(fc)] (A26) 

(A27) 



In the hterature, u; is used instead of Pq [ 4-6 ] . We may define the spectral functions 



and 



CTy (a;, T) = — Im Cy (a;, T) , 

TT 



ap{u!, T) — —Im Cp{u!, T) , 

TT 



(A28) 



(A29) 



Since different conventions are used in the hterature [4-6], we may use the notation 
ap{u!, T) and ay (a;, T) for the spectral functions given there. We have the following relations: 



ap{u,T)^ap{u,T), 



and 



av{uj,T) 3 



av{uJ,T) , 



(A30) 



(A31) 



2 4 

where the factor 3/4 arises because, in Refs. [4-6], there is a division by 4, while we have 
divided by 3, as in Eq. (A22). 
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APPENDIX B 



Here we extend the work of Appendix A to consider case of finite three-momentum, P. 
We consider the calculation of IiaJp{P^, P,T). The momenta P° and P are the values 
external to the loop diagram. Internal to the diagram, we have a quark of momentum 
k + P/2 leaving the left-hand vertex and an antiquark of momentum k — P/2 entering the 
left-hand vertex. It is useful to define 



E,{k) 



k + P/2 

(^k'^ + ^ + kP cos 9^ 



X 1/2 



and 



E2{k) = k-P/2 



P 



= k^ + — - kPcosO 



1/2 



Here k — \k\ and P 
We have 



Im Jy{P',P,T) = lN^(3ye{P') 



27r 



2Ei{k)2E2{k) 
x{[l - ni{k) - n2(A;)]5(P° - E^{k) - E^ik)) 

-[ni{k) - n2ik)]5iP' + Er{k) - E^ik)) 

-[n2{k) - ni(A;)]5(P° - E,{k) + E2{k)) 

-[1 - n,{k) - n2{k)]5{P'' + E^{k) + E2{k))) . 



Here, 



and 



ni{k) = 



n2{k) 



1 



1 



(Bl) 
(B2) 

(B3) 
(B4) 



(B5) 



(B6) 



(B7) 



e/3S2(fc) + 1 

In Eq. (B5), the second and third terms cancel and the fourth term does not contribute. It 
is useful to rewrite 6{P^ — Ei{k) — E2{k)) using 

2 



5[/(cos^)] 



I df I 
I d cos 6 I 



-5(cos^ — x) , 



(B8) 
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where 



We find 



= 


-- cos^e 














Of 


= \kp 




dcos9 


Ei{k)E2{k) 



(B9) 



and obtain 



Im Jp(pO,P,r) = \N,(5pe{P'){2^f j 



k^dk 



I 



2E^{k)E^{k) 



[I - ni{k) - n2{k)] 



(27r)3 
9/ (cos I 



(BIO) 



(Bll) 



9cos^ 
x5(cos^ — x)d{cos9) . 



We note there is a singularity when Ei{k) — E2{k). That occurs when cos^ = 0or^ = 7r/2. 
For our calculations we eliminate the point with 6 — t:/2 when evaluating the angular 
integral over d{cos9)5{cos9 — x) in the last expression. We obtain 



Im Jp(P°, P, T) = N^(5p e(P°)4^ / k'^dk e'^"!'^' 



(27r)3 

[l-r^i(A:)-n2(fc)] 



(B12) 



X 



kP\E^{k) - E2{k)\ 



where x is obtained from Eq. (B9), 



X 



kP 



p2 p2 

4 4 



1/2 



(B13) 
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